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ABSTRACT 


We  consider  the  perturbed  wave  equation: 


i  hi  at  rib-. t :  or,/ 
i  jv*) '  ,  *  ■; 


t.r.c.  or 


*  I  •  <  Jw  i 


utt  ”  +  e*(x#*#u)  *  o  for  x  e  ft  c  r 

t  e  r 

u(x,t)  “0  on  3(1  «  s 


where  £  depends  monotonically  on  u  and  is  periodic  in  t.  Periodic 
solutions  are  constructed  for  e  sufficiently  small  under  the  following 
hypotheses: 

(i)  (3  is  a  domain  such  that  the  eigenvalues  {X  }  of  -A  form  a 

n 

rational  sequence  -  i.e. 

\/*0  e  Q  . 

(ii)  The  period  of  f  is  rationally  related  to  the  periods  of  the  free 

2ir 

vibrations  for  the  wave  equation,  namely  - 

/T 

n 

(iii)  f  depends  monotonically  on  u  and  is  sufficiently  smooth  - 

fee  leads  to  a  solution  in  H  for  r  >  — r—  and  all  p  <  r. 


COfr 

1  "*S*>£CU0 

V  S  V 


This  generalizes  a  result  of  Rabinowitz  [6]  to  more  than  one  space  variable. 
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r~r? 


SIGNIFICANCE  AND  EXPLANATION 


\ 


r- 


In  [6] ,  Rabinowitz  proved  the  existence  of  2*-periodic  solutions  of  -the-  &- 


r 


one-dimensional  non-linear  wave  equation* — .'2—  y 


/ 


ufct  -  +  ef(x,t,u) 


(**) 


u(0,t)  ■  u(*,t)  “  0 

t  '  P‘- 

^for  sufficiently  small  and  f  2*^periodic  in  t ,  monotone  in  u,  and 
sufficiently  smooth.  This  answered  a  long-standing  open  question  and 
suggested  that  monotone  methods  could  be  used  to  overcome  solvability  problems 
in  bifurcation  situations  with  infinite  dimensional  kernels.  In  this  paper 
the  methods  of  Rabinowitz  161,  are  extended  to  higher  space  dimensions, 

r 

indicating  that  the  special  properties  of  the  one-dimensional  wave  equation 
are  not  essential  for  that  result.  What  remains  crucial  are  hypotheses  of 
rationality  in  the  relations  between  the  time  period  and  the  periods  of  the 
free  vibrations  for  the  wave  equation,  so  that  the  inverse  of  the  wave 
operator  remains  bounded  on  the  complement  of  the  null  space.  The  other 
crucial  factor  is  the  assumption  that  the  non-linearity  depends  monotonically 
on  u,  which  enables  us  to  solve  for  the  piece  of  the  solution  lying  in  the 
(possibly  infinite  dimensional)  null  space  of  the  wave  operator. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


PERIODIC  SOLUTIONS  OF  NON-DISSIPATIVELY  PERTURBED 
WAVE  EQUATIONS  IN  SEVERAL  SPACE  VARIABLES 


Robert  L.  Sachs 

1 .  Introduction 

We  consider  the  partial  differential  equation 
(*)  utt(x,t)  -  Au(x,t)  +  ef (x,t,u(x,t) )  -  0 

for  a  function  u(x,t)  where  x  is  an  n-vector,  t  is  a  scalar,  and  f  has 
period  2ir/w  in  t,  and  seek  solutions  satisfying  the  periodicity  and 
boundary  conditions: 


(1) 


u( x, t  +  2*/u)  -  u( x, t) 
u(x,t)  -  0  for  all  t 


whenever  x_j  ■  0  or  xj  “  Lj 


for  any  integer  1  <  j  <  n 


In  other  words,  the  x-variables  range  over  the  interior  of  a  rectangular 

parallelopiped  ft  e  Rn  and  we  impose  the  Dirichlet  boundary  condition  u  »  0 

on  3ft.  For  reasons  discussed  below,  we  restrict  our  attention  further  to 

those  ft  for  which  the  side  lengths  Lj  are  all  rational  multiples  of  one 

another  and  for  which  the  eigenvalues  of  -A  on  ft  (which  form  a  rationally 

2  2 

related  sequence)  are  rational  multiples  of  u>  .  Any  L  -function  v(x,t) 
satisfying  the  boundary  and  periodicity  conditions  ( 1 )  above  may  be 
represented  by  a  Fourier  series  of  the  form: 


(2) 


v ( x, t )  -  I  I 


ke*  jez" 


ikut 


sin 


’Vi 


sin 


if  j  x 
n  n 
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where  the  multi- index  j  e  z"  denotes  the  n- tuple  of  positive  integers 
( ji , j2»* # jn)*  If  Dv  ■  v  -  Av  is  in  ,  we  have: 


l  l 

kez  jezn 


-k  »  +  if 


4) 


ikoit  .  ^1*1 

a e  sin  — — ■■■  •  • . .  sin 
jk  I»1 


Our  rationality  restrictions  ensure  that  the  'Fourier  multiplier'  in  brackets 

£ 

choices  of  j  and  k,  is  bounded  away  from  0  for  all  other  choices  of  j 

and  k.  For  more  general  domains  in  several  space  variables,  -A  typically 

has  a  less  well-behaved  (in  the  sense  of  rational  multiples)  sequence  of 

eigenvalues  and  this  leads  to  number  -  theoretic  difficulties  of  'small 

divisors'  when  inverting  the  linear  wave  operator  □.  Even  in  one  space 

variable,  these  difficulties  have  not,  to  the  best  of  our  knowledge,  been 

overcome  to  date  on  problems  of  this  type.  In  the  case  where  -A  with  any 

given  boundary  conditions  has  such  a  spectrum,  our  results  hold. 

For  most  of  the  paper,  we  will  consider  n  =  2  and  Q  a  square  with 

side  length  «.  This  is  for  notational  and  expository  ease  and  the 

modifications  required  for  the  general  case  will  be  indicated  in  Section  5 

2  2 

below.  If  x,  y  denote  the  variables  in  the  square,  then  -A*  -3  -  9 
with  Dirichlet  boundary  conditions  has  eigenvalues  A  -  j  +  l  where  j,  l 
are  positive  integers.  Thus  2  is  the  lowest  eigenvalue  and  all  the  eigen¬ 
values  are  integers.  A  natural  choice  of  to  is  then  /i.  We  discuss  this 

case  and  again,  indicate  the  necessary  modifications  for  the  general  case  of 
2 

w  rational  at  the  end  of  the  paper. 


2  2  2 

above,  -k  <o  +  *  ( 


+...+  — ~),  although  it  may  vanish  for  infinitely  many 


Our  results  and  methods  are  the  direct  generalizations  of  Rabinowitz  [6, 
Part  I].  We  give  proofs  for  the  sake  of  completeness  but  acknowledge  our 
extreme  reliance  on  the  ideas  and  method  of  [6]. 

00  00 

Before  stating  our  results/  we  introduce  some  notation.  C  (C  will 

denote  the  infinitely  differentiable  functions  in  (x(y/t)/  ^2*  periodic  in 

t  (with  support  in  the  set  0  <  x,y  <  *).  Let  H®  *  L^  be  the  completion  of 
00 

C  with  respect  to 

I$Iq  ■  Ul2  -  I  ♦(x#y»t)  |2dxdydt 

a  oo 

and  H  be  the  completion  of  C  with  respect  to 

l$|2  ■  I  |D°$|2  where  o  =  (o  ,o  ,o  ) 

8  |o|«0  1  2  3 


>|o| 


lot  m  o^  +  °2  +  °3  wit*1  D  “ 


a.  o_  o, 

ax  1»y  2at 3 


and  let  R  denote  the  region 


0  <  x  <  w,  0<y<W/0<t<  /5».  We  denote  by  H8  the  completion  of  CQ 
with  respect  to  |  |g.  These  are  all  Hilbert  spaces  with  inner  products 
(  ,  )8.  Let  Cr  be  the  space  of  r  times  continuously  differentiable 
functions  of  x,y,t,  ^2  periodic  in  t  for  0  <  x,y  <  w.  we  denote  the  usual 
max  norm  in  Cr  by  I  -  thus 


l$l  *■  [  lD%l  where  l$l  “  sup  lt(x,y,t)|  . 

r  |o|-0  R 

Our  main  result  is  the  following: 

Theorem  1 .  If  f  e  d1  and  periodic  in  t,  with  f^  >  3  >  0,  then/  for 

|e|  sufficiently  small/  the  partial  differential  equation 
<*) 


u. .  -  u  -  u  +  £f (x,y,t,u)  “  0 
tt  xx  yy  * 


has  a  solution  u  e  hJ  such  that  u(x,y,t  +  /21O  =  u(x,y,t),  where 

2 

p  <  k.  In  particular ,  if  k  >  7/2,  u  is  C  and  hence  a  classical  solution 
of  (*). 

The  basic  idea  is  to  utilize  the  splitting  of  the  Hilbert  spaces  Hs 
into  two  pieces  -  the  part  intersecting  the  null  space  of  the  linear  wave 
operator  and  its  orthogonal  complement.  Here  the  null  space  is  infinite- 
dimensional,  so  this  may  be  viewed  as  a  Liapunov- Schmidt  procedure  in  a 
degenerate  case.  Thus  we  rewrite  (*)  as  the  pair  of  relations  with  u  **  v+ew 
where  v  e  N,  w  e  N^ 

f  n w  +  f (x,y,t,v+ew)  -  0  , 

(4)  / 

f(x,y,t,v+ew)  e  R(D)  -  N(O)  ^  . 

The  second  relation  is  the  'bifurcation  equation'  for  this  problem,  since  at 
e  »  0,  (*)  becomes 

(5)  Du  -  0  , 

which  has  infinitely  many  periodic  solutions.  As  in  [6] ,  the  bifurcation 
point  may  be  found  by  solving  a  convex  variational  problem,  which  we  do  in 
Section  2  below.  In  Section  3,  solutions  of  the  linear  wave  equation  are 
discussed.  We  note  that,  unlike  the  case  of  one  space  variable,  no 
derivatives  are  gained.  For  our  problem,  this  is  not  a  serious  difficulty. 
The  full  non-linear  problem  is  treated  in  Section  4,  while  in  Section  5  we 
discuss  the  modifications  needed  to  handle  more  general  periods,  unequal  side 
lengths,  and  higher  space  dimensions.  One  technical  argument  needed  in  the 
proof  of  Theorem  4  is  relegated  to  the  appendix. 


-4- 


2.  The  Bifurcation  Equation 


He  begin  by  considering  the  null  space  of  the  wave  operator,  □ ,  with 


the  boundary  and  periodicity  conditions  ( 1 ) •  It  is  easy  to  see  that  if  H  is 
the  L2-closure  of  the  smooth,  real  solutions  of  □  u  =  0  with  the  boundary  and 


periodicity  conditions  ( 1 )  - 


{u(0,y,t)  *  u(*,y,t)  *  u(x,0,t)  =  u(x,w,t)  -  0 
u(x,y,t  +  /2»)  ■  u(x,y,t) 

2 

then  N  =  {♦(x,y,t)  e  L  : 


♦(*#y,t)  -  l  ll  a.  te^iktsin  jx  sin  iy 

aW  3>°  £>° 


with  l  l  l  |aJja|2  <  -  and  *  g  - 


To  solve  (*),  split  u  into  pieces  v  e  N  and  ew  e  N^.  Then  u  * 

v(x,y,t,e)  +  ew(x,y,t,e)  with  Du  +  ef  ■  0  implies  that  we  must  have 
1  r 

f (x,y,t,v+ew)  e  N  .  Suppose  w  e  H  is  given.  He  shall  solve: 

(6)  f(x,y,t,v+Or)  e 

by  variational  methods  as  in  [6] « 

2  1 

Remark.  In  one  space  variable,  ♦  e  N  <*■>  $  ■  p(x+t)  -  p(-x+t)  p  e  L  (S  ) . 

For  several  space  variables,  $  is  not  so  simply  represented.  In  particular, 

2  2  2 

the  Diophantine  equation  2k  ■  j  +  i  has  many  solutions  besides  the 
obvious  sequence  j  *■  £  =  ±k  ■  s,  s  *  1,2,3,...  such  as  j  -  s,  4  »  7s, 
k  ■  ±5s  and  j  “  17s,  £  “  31s,  k  -  ±25s.  Nevertheless,  for  our  use  below,  we 
represent  ♦  e  N  as  follows: 


$(x,y,t)  "II  (b. j®1  3  *l  t  +  b  ,  e  4  3  +  t)sin  jx  sin  iy 
j>0  t> 0  3*  3* 

-  -  T  I  I  *  ♦  bTT  t)(ei(lw«y)  .  ,i(  jx-ty) 

j>0  t>0  3*  3* 

_  ei(-jx+iy)  +  ei(-jx-iy)] 


i  I  I  l*{b.t[ei(‘3xf|*W^  °  +  ei(jx"Ay+,/^  t]} 
j>0  A>0  3* 


■i  I  l  Re{b4#[ei(jXfty+/^  *  +  t,  }  -  ^(1)  _  ^(2)  # 

j>0  i>0  3* 


Note  that  each  term  extends  for  x,  y  beyond  the  square  and  in  this  extended 
sense,  each  3  ^  is  symmetric  about  the  origin  in  the  x,  y  plane;  also 
(-x,y,t)  ■  (x,y,t)  •  Thus  we  have 


M  -  max  (It1  |,  U  2  |)  ■  max  U1  | 
R  0<x<» 

-»<y<* 

0<t</2» 


0<x<w 

-»<y<W 

0  <t  <^2  w 


Also  observe  that  if  e  H*,  they  are  orthogonal  with  respect  to 

the  Hs  inner  product  for  every  s. 

He  shall  prove  the  following: 

Theorem  2.  If  f(x,y,t,u)  e  C2,  fu  >  0  >  0  and  w  is  given,  w  e  H2,  then 

there  exists  a  unique  v  e  N  n  H2  such  that 

///  f ( x,y,t,v+w)  ♦(x,y,t)dxdydt  -  0 


for  all  ♦  e  N  . 


Moreover, 


0/4  sup  |v(x,y,t) |  <  sup  |f(x,y,t,w)| 
R  R 


0|D  vL  <  C.  tlT  (x,y,t,v+w)  L  +  |f  (x,y,t,v+w)w^L] 


0|D2v|o  <  C2(|ftt<x,y,t,v+w)l0  +  2|fut(x,y,t#v+w)(vt+wt)l0 


where  ,  C2  are  constants* 

The  proof  of  Theorem  2  is  based  on  solving  an  approximate  problem  which 
is  better  behaved.  Namely,  consider  F(x,y,t,u) ,  where  Fu  “  f.  Let  Ng 
denote  the  subset  of  N  defined  by:  Ng  ■  {♦  6  N  |  |  $l2  <  K}.  Ng  is 
compactly  imbedded  in  N;  in  fact,  the  imbedding  is  a  Hilbert-Schmidt 
operator  (see  Adams  [1],  p.  174  or  Maurin  [3]).  He  shall  prove: 


Theorem  3. 


There  exists  Vg  *  Ng 


minimising  the  functional 


///  F(x,y ,t, f+w) dxdydt  over  all  +  e  NR  . 
R 


Moreover, 

j  P/4  sup  |vK(x,y,t)|  <  sup  |f(x,y,t,w(x,y,t)| 
1  R  R 


e(  l  |D°vk|2J/2<  C1Uft(x,y,t,VK4w>|()  +  |fu(x,y,t,vK+w)wtl0] 

|a|-1 

l  lDaVg|J)1/2<  C2[|ftt(x,y,t,vK+w)|0  +  2|fut(x,y,t,vK+w)(vt+wt)|( 


+  !0  *  '^U(  *'  *f  *,VK+W>Wtt 'o1 

t 


Proof  of  Theorem  3:  We  show  that  a  minimiser  vR  exists,  then  use  some 
particular  admissible  variations  about  Vg  to  derive  the  estimates  above. 

The  functional  /  F(x,y,t,  ^*w) dxdydt  is  continuous  on  Ng  in  the  H° 
topology  (w  is  fixed  and  bounded  and  ♦  e  NR  implies  ♦  is  bounded).  Hence 
by  the  compactness  of  Ng,  there  i.  .  sir  urn  which  is  taken  on  by  some 
Vg  e  Ng.  He  will  show  that  ,  when  K  is  sufficiently  large,  vK  is  an 


interior  minimum 


where 


To  this  end,  we  consider  variations  v„  +  n#  about  vv, 
v +  n#  e  N  for  h  <  0  sufficiently  small.  Then  we  have: 

&  K 

0  >  Jll  [F(x,y,t,v  +w)  -  F(x,y,t,v  +ri#H»)]dxdydt 
R  *  K 


=  ///  f (x,y,t,v  +w+np(x,y,t) #) (-n#)dxdydt 
R  K 


for  some  |p|  <  1  by  the  mean  value  theorem.  Dividing  by  -n  and  letting 
n  ♦  0,  np$  ♦  0  uniformly  the  'Euler  equation*  for  our  minimization  problem 
becomes: 

(7)  ///  f(x,y,t,v  +w) #  dxdydt  <  0  , 

R  K 

2  2 

for  all  admissible  variations  #.  If  |v  +n#|  <  K  ,  #  is  an  admissible 

variation  (vR+n#  e  Nk)  so  for  n  <  0  sufficiently  small,  it  suffices  to 
require  that  (VK'W2  >  °* 

0"e  admissible  variation  is  the  following: 


Let  q( X) 


0  if  |X|  <  m 

X-m  if  X  >  m 

X+m  if  A  <  -m  . 
..  „(1)  (2) 


Writing  VK  6  ^  as  ~  \  by  our  previous  remark,  choosing  m  »  M 


(vK,q(v”))  -  q(v«2,,,2  -  (v^1),q(v‘1)),2  ♦  <v<2,,(v<2)))2 
by  the  orthogonality  of  v*1  * ,  v*2*  in  H2.  Moreover, 


(v£i),q<v£i)))  -  l  (dV1’,  D^^v*1’))) 

K  K  2  ja|<2  K  K  0 


///  q(v''i))  +  q'(v<1))  l  lDVi)|2 

R  K  K  1  <|  a|  <2  K 


since  q" 


c  a.e 


From  our  choice  of  m,  v'  q(v'  ')  +  v'  q(v'  j  >  0  while  the  remaining 

K  K  Jv  K 

terms  are  non-negative.  Hence  q(vi.^)  -  q(v^,2S  is  an  admissible 

variation.  Substituting  in  (7),  we  have: 

///  f (x,y,t,vK+w) [q(v^ * )  -  q(v£2) )] dxdydt  <0  . 

R 

By  the  mean  value  thoerem,  we  have: 


///  [f (x,y,t,w)  +  fu(x,y,t, pvR+w)vR] [qtv^1  )  -  q(v£2  )]dxdydt  <0 


for  some  0  <  p  <  1 


Thus  we  obtain 


///  fu(int.  pt.Hv^1*  -  v^XqlVj*)  “  q(v^2)))dxdydt 

<  -  ///  f  (x,y,t#w)  [qfv*1  * )  -  q( yl2  ^ ) ]  dxdydt 

R  K  K 

<  sup  |f(x,y,t,w)  |  ///  (|q(v^1))|  +  |q(v^2 } ) | dxdydt  . 

R  R 

Using  the  monotonicity  of  q,  (vi.1  *  -  v^.2^ )  (qtvi.1  ^  )  -  q(vi.2*))  >0  so  that 
the  left-hand  side  is  bounded  above  (using  f^  >  3  >  0)  by: 

0  ///  )(q(v^1))  -  q(v^2) )) dxdydt 

=  0  ///  tv^qtv^1*)  +  v^2)q(v^2) )] dxdydt 

>  Bm  ///  ( I  q(v^  ^ )  |  +  |q(v*2) )  | )  dxdydt 
R  K  K 

since  Xq(  X)  >m  |q(A)| 


Thus 

Bm 

<  sup  |f (x,y,t,w) | . 

R 

But 

max 

|v  1  <  max  Ivl1 )  |  +  max 

(  2  ) 

|v^  |  <  2M  *  4m. 

D 

K'  _  K 

D  D 

Hence  max  |v  |  <  4m  <  4/3  sup  |f(x,y,t,w)| 


We  now  consider  a  different  variation.  Denote  by  cNx, y,t). 


the  time 


difference  quotient 


Ch(x,y,t)  =  £  (C(x,y,t+h)  -  C(x,y,t)] 


for  any  function  (.  We  will  show  that  the  centered  second  difference 
-  (v|j)  h  is  admissible. 

(V  "(VK)-h)2  “  (VK(t)'  ~2  (2VK(t)  ■  vK<t_h)  -  vR(t+h)))2 

h  v 

(rVR(t+h)  -  vR(t)  vR(t+h)  -  vR(t)^ 

“  ( - jj - «  - - )2  by  time  periodicity 

Ij  2 

*  lVK^2  *  ®  unless  vR  =  0,  in  which  case  the  subsequent 

estimates  are  trivial. 

Thus  by  (7)  again,  we  have: 

///  *(x*y#t,VR+w) l(-v£)  h]dxdydt  <  0  . 

R 

But  the  left-hand  side  is  equal  to  ('integrating  by  parts') 

///  <f(x,y,t,v  +w))h  v£  dxdydt 
R  an 


and  (f (x,y,t,vR+w) )h  -  ffc(int.  pt.)  +  fu(int.  pt.)(vK+w) 


h.  Thus 


///  f  (int.  pt . ) ( v  )  dxdydt 
R  u  K 


*  "  ///  (f..(int.  pt.)v£  +  f  (int.  pt.  )whvN dxdydt 
*  U  K 


which  implies: 


e,VKl0  <  Pt*,|o  +  lfu(int‘  Pt->whl0)  Iv^1 


K'O  * 


Letting  h  +  0,  we  have: 


0,VK  *0  <  ^0  +  sup  lfu{*'*'*'VK+w)|  *  lwtl0  • 

t  R 

But  ♦eN  implies  1^1*  -  •♦x*0+  *Vo’  1,1118 

B(,  _I  __  2(lft(,'''*'vK+w)l0  +  8UP  |fu(*#*,»,vK+w)|  |wt|Q) 


-10- 


Mow  consider  ♦  =  (((v^)  ^ ^ )  h.  By  a  calculation  similar  to  the 
previous  one, 

<v  ,W  ■  |-(vj))  h^2  >  0  (°r  else  t*ie  estimates  are  trivial)  . 

Substituting  ^  in  (7),  we  obtain: 


///  f(x,y,t,vR+w) [(v£)  h)h]  hdxdydt 
R 


///  _(tf(x,y,t,vK+w) {h)“h(-(v^)“h)dxdydt  <  0  . 


u  0  2 

Mow  (-f  )  is  a  difference  approximation  to  )  hence  the  mean 

value  theorem  implies 

(-fh)  h  -  -  ftt(»,*,«,int.  pt.)  -  2ffcu(int.  pt.)(v£  +  wh) 

-  fuu<int.  pt.)(v£  +  wh)2  +  fu(int.  pt.)(-(v£)”h  -  (wh)  h) 


8|“(V  *0  <  ///  Vint*  )2dxdydt 

<  ///  {[*tt(int.  pt.)  +  2ffcu(int.  pt.)(vR  +  w") 

+  f  (int.  pt.)(v!*  +  wh)2  +  f  (int.  pt.)(wh)  h)](-v^)  h}dxdydt  , 

UU  K  U  a 

whereupon  we  obtain 

0|-(v£)  hl0  <  |ftfc(int.  pt.)|Q  +  2jffcu(int.  pt.)(v£  +  wh)|Q 

+  |fuu(int.  pt.)<v£  +  wh)2|0  +  |fu(int.  pt.)(wh)“h|0  . 

Letting  h  +  0,  this  becomes: 


••Wo  <  lV','’'VV,l0  *  2|ftu<’'’','VW>(,Kt  *  Vo 


+  +  V  'o 
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Using  the  wave  equation  v  *  Av  ,  we  obtain  estimates  on  all  the  second 

order  derivatives.  This  completes  the  proo£  of  Theorem  3. 

Theorem  2  now  follows  rather  easily /  since  the  inequalities  of  Theorem  3 

imply  that  vR  is  an  interior  minimum  for  K  sufficiently  large.  In  this 

2 

case,  relation  (7)  becomes  an  equality  and  all  $  e  N  H  H  are  admissible. 

2 

Since  N  H  H  is  dense  in  N,  we  have: 

(9)  ///  f (x,y,t,v _+w) $  dxdydt  *  0 

R 

for  all  $  e  N.  Moreover,  vR  is  unique,  since  if  is  another  solution, 

0  -  ///  [f(x,y,t,vR+w)  -  f(x,y,t,vK+w)] (vR  -  vK>dxdydt 
R 

>  ei»K-  Vo  . 

which  is  a  contradiction*  Using  the  estimates  of  Theorem  3  we  obtain  Theorem 


Greater  regularity  of  v  can  be  obtained  by  assuming  more  on  f  and 
w.  In  particular,  suppose  f  e  Cr,  w  e  Hr  for  some  integer  r  >  2.  Then  we 
can  obtain  estimates  for  the  Hr-norm  of  v  and  some  pointwise  bounds  by  the 
same  methods  as  above.  For  instance,  using  the  variation  q(v  i1*)  -  q(v.(2)) 

V  t 

with  appropriate  choice  of  cut-off  point  m,  we  obtain: 

o  -  ///  [f  (x,y,t,v+w)  +  f  (x,y,t,v+w)(v +w  )] 

R  t  u  t  t 

[q(v^)  -  q(v£2* )]  dxdydt  . 

Thus  ...  .  . 

///  f  (x,y,t,v+w)v  (q(v^n ' )  -  q(v'2' )) dxdydt 
R  c  c 

<  ///  £-ft(x,y,t,v+w)  -  fu(x,y,t,v+w)wtJ [q(v^ } )  -  q(v*2) )] dxdydt 
which,  upon1  estimating  as  above,  leads  to: 


|  sup  |vfc|  <  sup  [|ft(x,y,t,v+w) |  +  |f  (x,y,t,v+w)w  |] 
R  R 
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Similarly,  using  higher  order  difference  quotients  in  time,  estimates  of  the 
form 


ll!t>P''lo  «  "J^o  +  TS  “Wo 

*  —  -  iv<-k>p»i„ ♦  i[,ls)P,1(vt4’t>lo 


which,  for  p  >  3,  using  the  Sobolev  inequality,  and  the  wave  equation 
readily  imply  that  v  e  Hr.  Indeed,  by  this  approach  one  proves: 

A  A 

Theorem  4.  For  2  <  p  <  r,  |v|  <  c  [|v|  .  +  |w|  +  1]  where  c  is  a 

r  r  f*  '  r  r 

constant  depending  on  f  and  its  first  p  derivative  and  depending 
monotonically  on  Iwl^  and  lv^_ I. 

For  a  full  proof,  see  the  appendix  below. 


3 


Inverting  the  Linear  Wave  Operator 


In  solving  the  nonlinear  wave  equation  (*),  we  will  need  to  solve  the 


linear  wave  equation: 


U  U  a  u 

tt  xx  yy 


f(x,y,t) 


IT  X 

for  given  f  6  H  ON,  with  w  satisfying  the  boundary  and  periodicity 
conditions  (1)  above.  The  solution  is  easy  using  Fourier  series. 

Theorem  5.  If  f  6Hrn  N'*',  31  w 

r  1  i 

□w  “  f  with  w  e  H  A  H  ft  N  . 

Proof.  Construct  a  weak  solution  in  H1  ft  N^  by  Fourier  series: 

f  “II  l  ®^iktsin  jx  sin  ty 

If  k  j>0  i>0  3 


2  2  2 
2k  -j  -4>0 


1  l  l 

k  j>0  i>0 

2  2  2 
2k-j-I>0 


/2ikt  .  ,  . 

e  sin  jx  sin  iy 


then  clearly  we  must  have 


**  -2kV+*2  * 

2  2  2 

The  term  -2k*  +  j  +  4  is  a  non-sero  integer  hence 
I (-2k2  +  j2  +  I2)  ^1  <1,  which  implies  that 

(wlj  <  | f | 1  and  more  generally 
|w|p  <  |f|p  for  1  <  p  <  r  . 

We  note  that  unlike  the  one  dimensional  case,  we  gain  no  derivatives  when 
solving  (10).  Indeed  the  sequence  j  ■  k-1,  l  -  k+1  yields 

-2k2  +  j2  +  42  -  -2k2  +  k2  -  2k+1  +  k2  +  2k+1  -  2 
thus  for  large  j,  k,  4,  y.^  decays  no  faster  than  at  least  for  this 

sequence  (and  other  obviously  related  ones). 

Note  that  here  our  rationality  conditions  seem  to  be  essential,  since  we 
have  then  good  control  on  the  norm  of  D-1  on  N^. 
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related  side  lengths,  the  spectrum  of 
'  problem  arises. 


4.  Solution  of  the  Non-linear  Problem 

We  wish  to  solve  (*)s  Du  +  ef(x,y,t,u)  *  0.  A  solution  of  the  form: 
u(x,y,t/£)  -  v(x,y,t; C)  +  ew(x,y,tj£) 

will  be  constructed,  where  v  e  N,  w  e  N^.  The  method  follows  that  of  [6]  and 
is  by  iteration* 

Proof  of  Theorem  1 :  Let  Ug(x,y,t)  =  u(x,y,t;e)  where  Uq  is  determined  by 
the  bifurcation  equation: 

f(x,y,t,uQ)  e  n\ 

Theorem  2  and  4  imply  that  for  f  e  Cr,  Uq  e  Hr. 

Let  w^  e  N'*’  solve  Dw^  "  -f (x,y,t,UQ ) .  Theorem  5  implies 
w1  6  N1  A  H* ,  and  then,  by  Theorem  2  and  4  again,  we  can  find  v^  e  N  O  Hr 
s.t.  ftx/y/t^+ewj)  e  S'*". 

Continuing  in  this  manner,  suppose  un_1  »  vn+1  +  ewjj_1  is  known.  Then 

|  y> 

solve  Dwn  -  -f(x,y,t,un-1 ),  obtaining  wr  e  N  OH,  and  then  find 
vn  e  N  O  Hr  from  the  bifurcation  equation  f(x,y,t,un)  e  NA.  If  {un} 
converges  to  u  in  a  strong  enough  sense,  we  will  have  constructed  a  solution 
to  (*). 


To  show  the  convergence  of  {u  },  we  follow  16)  and  first  assume  that 

n 

|unlr  is  bounded.  Define 


fiw  -  w  -  w  *  6v  -  v  -  v  ;  6u  =»  u  .«  -  u  *  6v  +  e6w 
n  n+i  n  n  n+ 1  n  n  n+1  n  n  n 

Then  Dfiw  -Ow  -  =»  -[f(x,y,t,u  )  -  f(x,y,t,u  , )). 

n  nr  in  n  n~  i 

Hence  |$wnlQ  <  |f(x,y,t,un>  -  f(x,y,t,u  )| 

*  I*.  <x,y,t,8(x,y,t)u  +  (1-0)u  )fiu  | 
u  n  n- 1  n-i  0 

<  sup  |f  (x,y,t,6u  +  (1-8)u  )|  |«u  |  . 

_  u  n  n- i  n- i  o 


Also,  f ( x,y , t ,u^} 


so  that 


///  [f(x#y,t,un)  -  £(x#y,t,un_1 )] 6vn-1dxdydt  -  0 

*  ///  f  (int.  pt.)[fiv  .  +  efiw  .]6v  .dxdydt  . 
u  e  n-1  n-1  n-1 

Thus 

e|6Vl'o  <  6  SUp  ,fuUnt*  Pt,)|  |fiwn-1!0 
R 


so  that 


|  6u  1 2  *  1 6v  |  2  +  e2 1  fiw  | 2 
nO  n  0  n  0 


<  e2(-^  SUP  lfuUnt.  pt.)|2  +  l)  l&fj2 


<  e2(—  sup  |  f  int.  pt.)|2  +  l)(sup  |fu(int.  pt.)  |  )2|  fiu^ 12 
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i.e.  Ifiunl0  <  6  (“J  B2  +  D^2B|  3uft_1l0  where  B  =  sup|  fu(x,y,t,p)  | .  If 

3 


|  e|  f— —  B2  +  1  )^2  B  <  1 ,  6u  ♦  0  geometrically  and  thus  u  u  in  H1 
a2  '  n  n 


By  interpolation ,  6u  ♦  0  in  HP  for  all  p  <  r  since 


Mu  I  <  C  |fiu  |J"p/r|6u  |P/r 
n  p  p  no  nr 


It  thus  only  remains  to  show  that  | r  is  bounded. 

For  this,  we  require  the  following  'composition  of  functions  inequality' 

[4]> 


(11) 


f  e  ck,  u  e  Hk  n  c°  , 

m 

|f(x,y,t,u(x,y,t)) \k  <  Ck(|u|k  +  1) 


where  C.  is  a  constant  depending  monotonically  on  lul  .  Indeed,  we  prove 

C° 

inductively  that  there  is  an  >  0  such  that,  for  some  positive  constants 


B,  K,  M, 


‘o'Vfl'r  <M'  ,vn'r  <K  ' 


If  Iv^  I  <  B,  Iv^j  lr  <  K,  eglwmlr  *  M  for  <  n  then  by  Theorem  2, 


Iv  I  <  4/0  sup  |f(x,y,t,ew  (x,y,t)|  . 

in  _  n 

R 

By  Sobolev's  inequality,  Iw  I  <  IwrI  <  alwnlr  *  <**  for  I el  *  eg*  Thus, 

C 

defining  x<«)  *  4/0  sup  |f(x,y,t,p)|  we  have 

x#y#t; |p|<s 

■V  «  *({■) )  • 

From  Theorem  4,  |vjr  <  cr t  e0  lwn,1 1  <  I vn lk_1  +  e0  lwnlk  +  11  for  I  el  <  e0* 

Since  |v  |  „  <  blv  l1y^r|v  |1-1/r  by  Sobolev  (for  r  large  enough)  we  have 
n  r- 1  n  i»  r 

A  W  I  . 

|v  |  <  c  [£„lw  ijfbf,  ”7  )1/r  |v  |  +  elw  |  +  l)  . 

n'r  r  0  n  1  v  v|v  |  J  n  r  0  n  r  1 

n  t 

By  considering  the  two  cases  where  the  factor  on  the  right  hand  side 
multiplying  |vnlr  is  <,  >  V2  respectively,  we  obtain 

|vn'r  <  2*rle0,Vl1(t0|,..lr  +  " 

A 

<  2crtOM) (M  +  1) 


| v  |  <  (2c  [ £_  Iw  IJb)rlv  I 

1  nr  r  0  n  1  n 

<  (2c  IoMjb)rlv  I  . 
r  n 

Adding,  we  have  in  all  cases: 

|v  |  <  2c  [cm] (M  +  1)  +  (2c  [oM]b)rlv  I  . 

n  r  r  r  n 

Finally,  from  Theorem  5  and  the  inequality  (11) 

|w  |  4  c  |F(x,y,t,v  +  EW  )  j  <  c  c  [  Iv  +6w  I]  ( |v  +0#  |  +  1) 

n+1  r  r  n  nr  rrnn  nnr 

<  C  C  [Iv  I  +  CM] ( | V  I  +  M  +  1)  . 

r  r  n  nr 

A  A 

Pick  ^(s)  >  max(X(s),  2c  [s] ,  (2bc  [s] )  ,  c  c  [s] )  monotone  non-decreasing 

r  r  r  r 

in  s.  We  have  derived  the  following  inequalities: 
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Iv  I  <  ♦(CM) 

n 

{  ,vJr  <  ♦<flM>(,vn  ■*•**+  1) 

l  '*W1 ‘  «V  *  H  ‘  "  * 

Pick  W*  1  (determing  M) ,  choose  B  -  1>(D  and  K  >  5(1)  (B  +  M  +  1)  * 
B(B  +  M  +  1 ). 

Then  the  first  two  inductive  steps  are  done  and  the  final  inequality 
becomes: 

|eraH.1  lr  <  *(B+1)»(K+M+1) 

so  choosing  eQ  «  M/1>(B+1 )  *{K+M+1 )  implies  c^lw^l^  *  **• 

To  complete  the  proof,  we  note  that  these  estimates  hold  to  begin  with: 
namely  wQ  “  0,  vQ  »  Uq  is  bounded  in  C°  norm  by  x(°)  since  w^  s  0  am 
X(0)  <  Xd)  <  t(1)  -  ».  Similarly, 

Jv0Jr  <  «1)(B  +  M  +  1)  -  K 

and  finally,  cq  <  CrCr *  ,V0  ^V0  ^r  +  ^  <  1KB)  (MH)  <  M  by  our  choice 

of  eQ.  This  completes  the  proof  of  Theorem  1. 

Mote  that  in  the  proof,  we  utilised  the  Sobolev  embedding  theorem.  Thui 
we  must  consider  large  enough  values  of  r,  namely  those  for  which 

Iw  I  .  <  a|w  |  for  some  constant  a  . 

n  c1  nr 

In  two  space  dimensions,  this  requires  r-1  >  3/2  and  for  N-space  variables 

N+1  taxi 

r-1  >  — .  if  r  is  slightly  larger  (r-2  >  — ) ,  then  since  our  solution 

P  2 

u  belongs  to  H  ,  p  <  r,  we  have  a  classical  (C  )  solution. 


5.  Modifications  for  Other  Domains#  Periods,  and  Dimensions 

The  differences  in  considering 
{*)  \»tt  -  Au  +  ef(x,t,u)  *  0 

with  boundary  and  periodicity  conditions  on  rectangular  parallelopipeds  with 

unequal,  but  rationally  related  sides  and  periods  rationally  related  to  the 

square  roots  of  the  sequence  of  eigenvalues  of  the  Dirichlet  Laplacian  are 

rather  minimal.  The  corresponding  divisor  when  solving  Dw  =  g  for  w  e  W*- 

is  a  rational  number  bounded  away  from  zero,  hence  □“*  is  well-defined  and 

bounded  (perhaps  by  a  constant  larger  than  1).  Similarly,  when  using 

embedding  theorems  for  several  space  variables,  the  assumed  smoothness  of  f 

must  be  correspondingly  higher.  The  estimates  of  Section  2  are  identical, 

except  that  Nj^  must  be  a  closed  ball  in  the  H  P-norm  for  p  >  N+1/2  and 
a 

f  must  be  C  •  The  same  trial  variations  of  vK  will  work,  although  more  t- 
derivatives  must  be  used  if  N  >  4. 

The  result  also  extends  to  other  domains  and/or  boundary  conditions 
provided  that  the  eigenvalues  of  -A  satisfy  the  required  rationality 
conditions.  Thus,  for  instance,  we  could  consider  rectangular  parallelopipeds 
with  rational  sides  and  impose  periodicity  in  each  x-coordinate  separately  or 
a  combination  of  periodicity  and  Dirichlet  boundary  conditions.  Of  course, 
for  several  dimensions,  the  rationality  restriction  is  more  severe  than  in  one 
space  variable.  What  remains  crucial  is  that  -  XJ  >  v  >  0  for  all 

i  >  0  where  X^  are  the  eigenvalues  of  -A,  excluding  those  values  l,  n 
for  which  this  vanishes  (the  null  space  of  ). 

We  note  also  that  the  splitting  of  N  into  two  orthogonal  pieces  - 
♦  e  N  <*■»>  $  ■  ^  which  we  used  in  making  several  key  estimates, 

can  be  accomplished  for  any  number  of  space  variables  (for  an  odd  number  of 
space  dimensions,  the  expressions  involve  the  imaginary  parts  of  sums  cf 
exponentials  instead  of  real  parts,  but  are  otherwise  the  same). 


Appendix 

Proof  of  Theorem  4:  (See  [4],  [6].)  Suppose  v  e  Hr  A  N  with 

f  (x,v+w)  e  N"^  for  given  w  e  Hr  (here  x  e  0  c  Rn  and  the  time  period  is 

T).  We  wish  to  show  that  |vj  <  c  [Iwl.,,  lv IJ(|v|  ,  +  |w.|  +  1). 

r  r  i  t  r- 1  r 

Since  v  e  N,  it  suffices  to  estimate 

|3*v|  <  crt»wl1,  lvtl](|v|rH  +  |w|r  +  1) 
for  then  we  have  control  over  all  r**  order  derivatives,  and  lower  order  terms 
will  be  readily  absorbed.  Now  f  e  s\  u«v  +  weH_AC  implies 
3*f (x,y,t,v+w)  e  N^.  Let  ffca  denote  a  derivatives  of  f  with  respect  to 
the  third  argument  t,  and  similarly  for  fu«  Since  3*v  e  N,  we  have: 

0  *  //  3*f (x,y#t,v+w)  •  3*v  dxdydt 

Qx[0,t] 

■  //  t(f  (x,y,t,u)  (3*v  +  3*w)  3*v]  dxdydt 
nx[o,T]  u  t  t  t 

+  //  [3*f  -  f  (3%  +  3%)]  (3%) dxdydt  . 

a*tO,T] 

Hence,  by  monotonicity  and  f^  >  B  >  0,  we  have: 

•  l<v|0  <  lfu<3tv  +  3tw)  ‘  3tf|0 

+  sup  |f  (x,y,t,p)|  •  I  3  w | 

x,y,ten,t+0,T] 

|p|  <L 

where  L  is  the  bound  on  sup  u  from  the  Sobolev  inequality.  Thus  to 

complete  the  proof,  we  show,  as  in  [6] ,  that 

1 3fcf  (x,t,u)  -  fu  3'u|q  <  crtlwl1,  «vtl](|v|r_1  +  ♦  1) 

which  holds  in  general  for  Cr  functions  f.  The  difference  consists  of 

a  .  ol  a 

terms  ..  the  form  c  f  ( 3.  u)  1  ( 3“u)  2  •••  Of”  u)  r~  with 

°Pa  t  t  t 
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+  a2  +•••+  o^  -  p  and  r  -  o  -  +  2a2  +•••+  (r-Da^  and  c  a  »« 

r_1  i 

positive  binomial  type  coefficients.  Consider  f  n  ( 3d)  -  let  Uq  * 

t  uP  3-1 

f  and  u*  *  3^u,  j  -  1,...,r-1.  For  j  >  1#  we  choose  p.  -  77^77 ■  . 
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.  To  use  the  Hfllder 


inequality  on  the  above  product,  we  wish  to  choose  ~  -  1  -  (“^2)  -  ~f~~2 

which  is  fine  for  0  +  p  >  2.  For  p  -  0,  the  given  term  is  f  _,  which  we 

tr 

bound  by  a  sup  norm  since  u  is  bounded  by  assumption,  p  -  1  leads  to 

a-0 

JT 

the  only  term  we  exclude  in  our  sum,  namely  f  3  u.  Thus  all  other  terms 

u  t* 

satisfy  0  +  p  >  2  and  we  choose  pfl  as  above.  Then  Holder's  inequality 
gives: 

r-1  2o.  r^1  2op  VPa 

//  Hu3  dxdydt  <  H  (//  u  3  3  dxdydt)  3  . 

flx[0,T]  j-0  3  j-0  flx[0,T]  3 

Applying  the  Gagliardo-Hirenberg  inequality  to  a  typical  term  in  the  product. 


iff  (33u)2(r'2,/(j-1)  dxdydt]  <J-3>/2<-2)  .  .A,  „ 
ax[0,T]  t  L  1 


r-2)/3-1 


.  .  1- ( j-1/r-2 ) ,  , j-1/r-2  .  . 

<  const.  Iufcl0  |u|J_1  +  |utlQ  . 


But  I u^t 0  <  (vo*  R)'2  lut*0  and  |ufcl0  <  •ulr.1»  80 

US  (a3o)2(r-2)/(j-1))(j-1)/2(r-2)  ,  c,u  ,1-< j-1 >/(r-2>  (j-,)/(r-2) 

flilO.T)  *  1  0  r'1 


Thus  the  product  to  be  estimated 


1/p.  r-1 


3-0  ox[o,t]  3  j-i 


which  proves  Theorem  4 
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